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Abstract 

The derivative of self-intersection local time (DSLT) for Brownian motion was in- 
troduced by Rosen HRos05] and subsequently used by others to study the L 2 and L 3 
moduli of continuity of Brownian local time. A version of the DSLT for fractional 
Brownian motion (fBm) was introduced in HYYL081 . In the present work, we further 
its study by proving a Tanaka-style formula for it. In the course of this endeavor we 
provide a Fubini theorem for integrals with respect to fBm. The Fubini theorem may 
be of independent interest, as it generalizes (to Hida distributions) similar results pre- 
viously seen in the literature. Finally, we also give an explicit Wiener chaos expansion 
for the DSLT of fBm. 
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1 Introduction 

Let B t denote a Brownian motion on R with Bq = and let L(t, x) = Lf be its local 
time at x up to time t. In connection with stochastic area integrals with respect to 
local time [Wal83] and the Brownian excursion filtration, Rogers and Walsh |R W91al 
studied the space integral of local time, 

ft QJ± 

A(t,x):= J l [0:Oo) (x-B s )ds ; — = L(t,x). (1.1) 

In a companion work on the Brownian local time sheet [RW91b], the occupation 
density of A(t,B t ) was investigated. It was also shown in BRW901 that the process 
A(t,B t ) — J Q L(s, B s )dB s has finite, non-zero 4/3-variation. An alternate proof of 
this fact using fractional martingales was recently given in [HNS12]. 

In [Ros05], Rosen developed a new approach to the study of A(t, B t ) as follows. 
If one lets h{x) := \o,oc){ x )> tnen formally 



d d 2 
—h{x) = Six) and -r^hix) = S'(x), (1.2) 
dx dx z 



where 5 is the Dirac delta distribution. Holding r fixed and applying Ito's formula 
with respect to the Brownian motion B s — B r gives 

l [0jO o)(^ - B r ) - l [0 ,oo)(0) = J 5{B S - B r ) dB s + I J 5'(B S - B r ) ds. (1.3) 

Integrating with respect to s from to t and interchanging the orders of integration in 
the integrals leads to 

A{t,B t )- [ L{s,B s )dB s = t+\ f [ 5'{B s -B T )drds. (1.4) 
J ^ J J 

This formal identity was stated in [Ros05|. We note however, that there is some am- 
biguity since if we change the definition of h only slightly to h(x) := l(o,oo)( x ) an d 
apply Ito's formula in the same manner, we obtain 

A(t,B t )- [ L{s,B s )dB s = \ I I 5'{B S -B r )drds. (1.5) 
JO 1 Jo Jo 

Motivated by (11.41) . Rosen |[Ros05ll showed the existence of a process now known 
as the derivative of self-intersection local time (DSLT) for B t . That is, he demonstrated 
the existence of a process a' t (y) formally defined as 

<(y)-=- I I 5'{B s -B r -y)drds. (1.6) 
Jo Jo 
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This process was later used in [HNQ91 and HHN101 to prove Central Limit Theorems 
for the L 2 and L 3 moduli of continuity of Brownian local time. 
In [Mar08a], it was proved that almost surely, for all y and t, 

\ oc'M + ~ sgn(y)t = J Lf°- y dB s-\J^ sgn(B t - B r -y)dr. (1.7) 

In particular, for y = equation (11.51 ) is correct and (11.41 ) should not have a t term. 
The formula (11.71 ) is commonly referred to as a Tanaka formula. The method of proof 
in MMar08al also serves to give an alternate proof of the existence of a' t {y) and joint 
continuity of a' t (y) + sgn(y)i, which Rosen had deduced earlier by other methods. 
In MMarl2ll . yet another existence proof for a' t (y) was given using its Wiener chaos 
expansion. Our aim is to extend these results to the more general process of fractional 
Brownian motion. 

Standard fractional Brownian motion (fBm), with Hurst parameter H 6 (0, 1), is 
the unique centered Gaussian process with covariance function 

E(B?B?) = ±(s 2H + t 2H -\t- s\ 2H ). (1.8) 

Note that H = 1/2 gives us a standard Brownian motion. In HHuOU . it was shown that 
the self-intersection local time of fBm is differentiable in the Meyer- Watanabe sense. 
Using Hu's arguments, for H < 2/3, [YYL08] deduced the existence of processes 
which are related to what we shall call the DSLT of fBm. We should note that a 
slightly erroneous bound is utilized in both [HuOH and HYYL081 . for which we have 
provided a corrected modification in Appendix |A] 

The work of [YYL08] shows that there is a critical value below which the p- 
variation of their DSLT of fBm is non-trivial. Their result was partially motivated 
by stochastic integrals with respect to fBm-local times which were further studied in 
MYLY0911 . In particular, in HYYL081 . two versions of a DSLT of fBm are defined and 
the following formal identity is stated 

Ha' t (0)=t+ f Lf" dB* - \ f sgn(£f - B?) dr. (1.9) 
Jo z Jo 

This reduces to ([L4l when H = 1/2. 

In this work, by proving an analog of (11.71 (which again shows that (11.91 ) is better 
off without a t term), we modify equation (11.91 to a process formally defined by 

a' t (y) := -H f f ' 5'(B? - B? - y){s - r) 2H - 1 drds. (1.10) 
Jo Jo 

When H < 2/3, we will show that such a process exists in the L 2 sense. 
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The rest of the paper is organized as follows. In the next section we give a few 
remarks on the existence of a' t (0). In Section [3] we review some tools from Malliavin 
calculus needed for the sequel. One of these tools is a Fubini theorem for fractional 
Brownian integrals which generalizes, to Hida distributions, similar results found in 
[CN03 and IIMis081 . In Section 0] we present an explicit Wiener chaos expansion for 
a' t (0) and prove the existence of DSLT for all y G R. We conclude in Section [5] by 
proving a Tanaka formula for a' t (y). 



2 Existence of the DSLT of fBm 

Let FoG denote the Wick produclQ of F and G. Formally applying Ito's Lemma (see 
Theorem 13 . 8 I below) for fBm in a similar fashion to (11.51 gives 



l (0lOo) (Bf - Bf) = J* 6(B? - B?)«dB? + hJ* S'(B* - B*){s - r) 2 ^ 1 cfa. 
Integrating with respect to r and switching order of integration gives 

[\ i0tOo) (B?-B?)dT= f L^odBf+H f I* S'iBf-B^s-rfX-Urds. 
Jo Jo Jo Jo 

(2.1) 

Here (s — r) 211 ^ 1 is called a reproducing kernel, and is standard in integrals involving 
Bf (see Chapter 2 of ||BH0ZO8B ). Comparing the above with (fL4l and (fL6l >. it is 
natural to define the DSLT of fBm as the formal process a' t (y) given in (1 1 - 10b . 

In order to rigorously define a' t (y), let /i( 2;) denote the standard Gaussian density. 
Weset/ £ (x) := ^f(f e ), so that * 

Mx) = -^e-^ 2 /e and f' e [x) = ±f 6 {x). (2.2) 

Note that as e — > 0, f £ (x) converges weakly to the Dirac delta distribution, 5(x). For 
fixed < H < 1, let 

<x' t ,M := ~ f f ' f' e {Bf -B? -y)( s -rf H ~ l drds. (2.3) 
Jo Jo 

An analog of the following result appears in HYYL Q8 1. 

Proposition 2.1. For fBm with H < 2/3, defined on the probability space (Q,J-,V), 
the processes a' t e (0) converge in L 2 (V) to a process a' t (0) as e — > 0. Moreover, a' t (0) 
is continuous in t. 

'The Wick product is further explained in Section [5] or for a full treatment see [BH0ZO8 1 . 
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The proof of the above is similar to the arguments of HYYL081 . except for the 
convergence of an integral which is given in Lemma l2~2l below. The proof presents 
some points of interest, so we sketch it as follows. 

The key lies in computing E[(a' t e (0)) 2 ]. In the sequel, K denotes a positive con- 
stant which may change from line to line. We start by expressing f £ (x) in a convenient 
form using the Fourier identity f\{x) = -h- L e ipx e~ p I 2 dp. This gives 



fe(x) 



1 



whence 

Now, for V t = {0 < r < s < t}, 



? ipx e- £p2/2 dp, 



i 

2^ 



pe 



ipx p -ep 2 /2 ^ 



(2.4) 



(2.5) 



K 



pqe 



-e(p2W)/2 p -Vai(p(B?-B?)+ q (B^~B%))/2 



T>? J(p,q 



x s — r 



,2H-\t I 



(s' - r'f u - x dp dq dr dr' ds ds' . 



We will show that this can be bounded uniformly in e. Using standard notation from 
the literature, let 



A 


:= Var(Bf-5f) 


= \s-r\ 2H 


P 


:= Vai{B$-B$) 


= \s'-r'\ 2H 


/' 


:= Cov(#f - By 


B s r — B r i ) 




= \(\s'-r? H + 


s — r \ — , 



s' - s\ 2H 



\r'-r\ 2H ). 



With this notation, we have 



E[« £ (0)) 2 ] = K [ {s-r? H -\s>-rT- 1 
Jv 2 



pqe -e(p 2 W)/2 e -(p 2 X+2 Pq „W P )/2 dp dq dr dr > ds ^ 



Isolating the dq integral we have: 



qe -pq^ e -Q 2 (p+e)/2 dg 



e 2(p+e) 

2 2 
P M 
g2(p+e) 



-(q+ J ^) 2 (p+e)/2 , 

qe p+ € ' " dq 



qe 



-q 2 (p+e)/2 



dq 



p/j, 



(p + e) 



-q 2 (p+e)/2 _ 



(2.6) 



(2.7) 



(2.8) 
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The first term on the right side is an integral of an odd function, and thus vanishes. 
The second integral on the right side converges as e — > to 

Kr V 2 ^/(2p) W 

At e = 0, we therefore have for the dp integral, 

Hjt f 2 e ^(Ap-M 2 ) d = 

pW/ P (Ap-M 2 ) 3 / 2 ' 

Thus, we have reduced the problem to determining the integrability, over V 2 , of 

K^{s-r) 2H - l {s' -r') 211 - 1 



(Xp - p?fl 2 



(2.9) 



The existence of a' t (0) is therefore proved by invoking the following lemma, which is 
proved in the appendix. 

Lemma 2.2. If H < 2/3, then 

r a ( s _ r )2H-l( „/ _ r ')2H-l 

/ drdr' dsds' <oo. (2.10) 

j<d\ (A/? - p 2 yi 2 

A few remarks are in order. First, something close to Lemma 12.21 was proved in 
MYYL 081. but they had a factor of s 2H ~ l where we have (s — r) 2H ~ l . This resulted 
from applying Ito's formula to Bf as opposed to Bf — Bf. 

Next we note that in [Ros05], Rosen states "The (DSLT of Brownian motion) in 
R 1 , in a certain sense, is even more singular than self-intersection local time in R 2 ." 
If we believe that the critical Hurst parameter H c for the DSLT to exist in L 2 is 2/3, 
then Rosen's statement would be supported by the fact that 2/3 is less than the critical 
Hurst parameter for the self-intersection local time of planar fBm (H c = 3/4, see 
nRos87llHN031 ). 

Above the critical parameter H c , the behavior of a' t e (0) as e — > is also of in- 
terest. One would expect a Central Limit Theorem to exist, along the lines of Theorem 
2 in HHN051 or Theorem 1 in [Mar08b ], but this remains unproved. In particular, it 
seems as though the techniques developed in [IHN05M should apply, especially since the 
Wiener chaos expansion for a' £ is readily computed (see Section @]), but the presence 
of the derivative seems to complicate matters. Nevertheless, we venture the following 
conjecture. 

Conjecture: 

• The critical parameter is H c = 2/3. At H c , log (jy e ) 7 °4, £ (0) converges in distri- 
bution to a normal law for some 7 > 0. 
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• For H > H c , £-^a{ £ (0) converges in distribution to a normal law for some 
function j(H) > which is linear in 1/H and for which 7(2/3) = 0. 

This would mirror the behavior of the intersection local time as seen in IIHN051 . We 
should mention that, to our knowledge, no such Central Limit Theorem has yet been 
proved even for intersection local time in two dimensions at H c = 3/4. 

3 A Fubini theorem for the WIS integral 

There are several different ways one can integrate with respect to fBm as can be seen in 
[BH0ZO8] or [Mis08]. We use the integral based on the fractional white noise theory 
introduced in HEVDH03H (s ee also HB0SWO4H , HBH0ZO81 Ch. 4]). In particular, we 
adopt the nomenclature of [BH0ZO8] and call this stochastic integral the Wick-Ito- 
Skorohod (WIS) integral. 

In an effort to be somewhat self-contained, in this section we summarize some re- 
sults concerning white noise and the WIS integral. For more details we refer the reader 
to HHQUZ96I IEVDH031 |B0S W04B . The only new result in the section is Theorem ET9l 
which is standard for other stochastic integrals (cf. [CN05, Theorem 3.7], |Mis08 ( 
Theorem 1.13.1]); however, we have found no reference for such results with respect 
to Hida distributios and the WIS integral. Theorem 13.91 follows easily once one has 
the right definitions. Its formulation may prove useful in its own right, but the main 
reason for its presentation here is due to its role in proving the Tanaka formula. 

3.1 Classical white noise theory 

Let A = Nq be the set of multi-indices with finite support, <S(R) be the Schwartz 
space, and S'(R) be the space of tempered distributions. By the Bochner-Minlos 
theorem there is a probability measure V on the Borel (T-field of S' satisfying 

J s ^e X p(i(uJ))dP(u)=e^(-^\\f\\l^ , fES(R). (3.1) 

This measure satisfies, for all / G <S(R), 

E(oj, f) = and E(w, f) 2 = ||/||£ 2(K) . (3.2) 

If g n S <S(R) converge to l [0)t] in L 2 (R), we define (u, l [0)t ]) := lim^^u;, g n ) as a 
limit in L 2 {V). The family {^[o,t]) t>0 now maps <S'(R) to R, and can be identified with 
random variables having characteristic functions (j){v) = exp (— i(ti^) 2 ). Choosing a 
continuous version of this family gives us Brownian motion. 
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Recall that the Hermite polynomials given by 

K{x) := (-lfe* ' 2 — (e~ x I 2 ), n = 0, 1, 2, . . . (3.3) 

are orthogonal with respect to the standard Gaussian measure on R. Thus multiplying 
by the Gaussian density and choosing a convenient normalization gives us an orthonor- 
mal basis for L 2 (R), 

e n (x) := vr- 1 /^ - l)\)- l l 2 h n ^{V2x)e- x2 l 2 , n = 1, 2, 3, . . . , (3.4) 

called the //enra'te/wTicf/oft.s'H 

For /3 = . . . , j3 n ) £ A, we define 

^(w) := ^((a;,ei})%({a;,6}) ■ ■ ■ h^u, £ n )). (3.5) 
Every G L 2 {V) has a representation in terms of the 7^: 

F( W ) = 5]c^( W ) (3.6) 

/3GA 

where the series converges in L 2 (V). Moreover, one has the isometry 

bf 2 = y,^4- 

/3eA 

Representation (13.61 ) for F is called the Hermite chaos expansion, and it is related 
to the Wiener-Ito chaos expansion^] via the following formula which follows from 
MTU : 

f e P dB &n = Up{u>). (3.7) 

Here denotes the symmetrized tensor product, and ^ 0/3 := • • • ■ In 
particular, the Hermite chaos is a way of writing the nth Wiener-Ito chaos in terms of 
n-fold products of Hermite functions, which form orthonormal bases of L 2 (R ra ) (see 
HHOUZ96I pg. 30]). 

The main reason for using chaos expansions in terms of Hermite polynomials in- 
stead of multiple Wiener-Ito integrals is the natural extension to distributions they have 
from L 2 (V) random variables. Let 

(2N)T:=(2-ip(2.2)'»...(2-nP» for 7 = ( 7l , . . . , 7n ) e A. (3.8) 

2 One may substitute any orthonormal basis of L 2 (R) whose elements possess decay properties such that 
Lemma 4. 1 of EVDH031 holds. See also Theorem 3. 1 in lHto3Tl . 

3 For a full treatment of multiple Wiener-Ito integrals and their related chaos expansions, see [Nua95 1. 
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Definition 3.1 (Hida test functions and distributions). Given the probability measure 
V on S 1 , the Hida test function space (S) is the set of all tp G L 2 (V) given by 

/3eA 

satisfying 

J^4/3!(2N) fc/3 for all k = 1,2,... 

p<=A 

The Hida distribution space (S)* is the set of all formal expansions 

/3eA 

satisfying 

^/3!(2N)~ 9/3 for some q€R. 

/3eA 

It was shown in RZha921 that (S)* is the dual of (S). Moreover, by Corollary 2.3.8 
of HHOUZ96L 

(5) C L 2 (V) C (S)*. 

This should be thought of as analogous to the triplet S C L 2 (R) C S*. Note that if 
F e L 2 (P),then 

((F^)) = (F^) LHP) =B(FiP). 

Thus, for ip = Yy/3eA a ^l 3 e ("->) anc * ^ = S,aeA g (^)*> tne duality inher- 

ited from L 2 (V) is given by 

((F,^)):=J2^h- (3-9) 

/3eA 

Lemma 3.2. Suppose F{x) is an (S)* -valued function on a a-finite measure space 
(X, £>, v) and that 

((F(x),4>)} £ L\X,v) VVG(<S). (3.10) 
77je« f/zere a unique G in (S)* such that 

«G,V))= / {{F{x),ip))dv V^e(5). (3.11) 

We write J x F(x) du := G. 

Proo/ See Theorem 3.7.1 in IIHP571 or Proposition 8.1 in HHKPS93L □ 
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Let T C R be a time interval. In light of the above result, we say an (5)*-valued 
process is in L l (T) if 

G L\T) for all V € (S). (3.12) 

The Hida distribution space (<S)* is a convenient space on which to define the Wick 
product. 

Definition 3.3 (Wick product). If F = J2/3eA bpUp and G = 2~2^eA c 7%/3 are de- 
ments of (S)*, then their Wick product is defined as 

FoG := } ] bpCyHp+^ju) 

13,-7 

= E ( E bpcJnxiu). 

AeA \/3+7=A / 
By Lemma 2.4.4 in [HOUZ96], (S)* is closed under this product. 

Given the measure V on S' from (13. It . we may write Brownian motion as 

B(t) = (w,l [ o, t] (0) = /w,^(l [ o, t ],e fc ) ia( R ) efc(o\ (3-13) 
= E / £,k{s)ds(u,£, k ) = E/ 6(s)(is'H e (fc)(w). 



fe=i ,,u fe=i' 



where eC ) is the multi-index with a 1 in the fcfh entry and O's elsewhere. This motivates 
the following definition: 

Definition 3.4 (White noise). The (S)* -valued process 

W(t) :=£ek(iRw(w) 
k>i 

is called white noise. 

Using the Wick product on (S)* and the definition of the integral of a Hida distri- 
bution given in Lemma [3721 we can now integrate a Hida distribution with respect to 
white noise: 

/ F(t)dB t := [ F(t)oW(t)dt. (3.14) 
Jr Jr 

The above is called a WIS integral with respect to white noise. The following theorem 
shows that the WIS integral is a generalization of the Skorohod integral: 
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Theorem 3.5. Suppose F(t,cj) : Rxfi — >• R is Skorohod integrable. Then F(t, •) o W(t) 
is dt-integrable in (S)* and 

rb rb 

/ F(t,u;)SB(t) = / F(t,-)oiy(t)(it. (3.15) 

Ja Jo 

Proof. See Theorem 2.5.9 in IHOUZ961 . □ 



3.2 Fractional white noise theory and a Fubini theorem 

Elliot and Van Der Hoek [EVDH03] introduced the fractional white noise as an ele- 
ment of the Hida distribution space, and thus constructed the WIS integral which is 
valid for any H G (0,1). The main tool used to define the fractional white noise is the 
following operator for which we set: 

c H ■= v / sm(7ri?)r(2F + 1) 
Definition 3.6 (Mjj operator). The Mh operator on f G <S(R) is defined by 

M^f(y) = c H \y\ 1 ~ H f(y). 

This operator extends to the space 

L 2 H (R) := {f:M H feL 2 (R)} (3.16) 
= {f:\y\^ H f(y)eL 2 (R)} 

which is equipped with the inner product 

(/^)l^(r) = (M H f,M H g) L 2 m . (3.17) 

We note that L 2 H (R) is not closed under this inner product, and that its closure contains 
distributions (see HNua95l pg. 280] or HBH0ZO81 Ch. 2]). 

Since M H f £ L 2 (R) we can moreover define M H : <S'(R) -)• <S'(R) by 

(M H u,f):=(u,M H f) for/G5(R),^e5 , (R). (3.18) 

We have defined Mh for a given fixed H G (0, 1), however, the theory extends to M 
operators which are linear combinations a^M^ + • • • + a n M^ n ; for more on the M 
operator see ISKM871lEVDH551ll^lTI . 

4 Compare this to the integral of [H0O3] where a fractional Brownian measure is defined directly on 
tempered distributions S' . Integrals with respect to this measure are defined for H > 1/2. 
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Using the orthonormal basis e& := M H l ^k of L 2 H (R), we define (and choose a 
continuous version) of fBm as 

B H (t) := (oj,M H l m (-)) = (M H oj,l m (-)) (3.19) 

which satisfies, by (O and (A. 10) in HEVDH03L 

E[B H (t)B H (s)] = (M H l m ,M H l [0>s] ) Lm (3.20) 

= l {t 2H + s 2H_ lt _ sl 2Hy 

Note that the action of S' on S given by (ui, •) is still inherited from L 2 (R) and not 
from L 2 H (R). 

The definition in ( 13. 191 ) can be further rewritten as 

B H (t) = ={M H u,l M (-)) (3.21) 
= ^Mhw, J^(l[o,t],e fc ) i ^ (R je fc (-)\ 

/ oo \ 



^ M ff w , i M H 1 [0,t] , ik > L2 (R) efc ( • ) J 

OO CO 

V/ M H ik{s)ds{M H uj,e k ) = V/ M H £ k (s)dsH t 
k=i Jo k=i Jo 



(fc) (w) 



which motivates the following notion of fractional white noise: 

W H (t) := Y M H ik(t)U e(k) (w). (3.22) 



fc>i 



By Lemma 4.1 in [EVDH03], W-^i) is an (5)*-valued process. We note here that 
the underlying probability measure V on S' is the same as for W(t). 

Definition 3.7 (WIS integral). Let F(t) be an (S)* -valued process such that F(t) o W H ( 
jj iw L 1 (M) (as in (13.121 ) ). We define 

[ F{t)dB? := / F{t)odB?:= [ F(t) oW H (t) dt. 
Jr Jr Jr 

Theorem 3.8 (fractional Ito formula). Let f(s, x) : R x R ^ R be in C 1,2 (R x R). 
7f random variables 

f(t,Bf), f f s (s,B?)ds, and f* f xx (s, B?^ 21 *' 1 ds 
Jo Jo 
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are in L 2 (V) then 

f(t,B?)-f '(0,0) (3.23) 

= I f B ( 8 ,B?)da+ [ f x (a,B?)dB? + H [ f xx {s,B»)s 2H ~ l ds. 
Jo Jo Jo 

Proof. See Theorem 3.8 of HB0SWO41 . □ 

We now prove a result which one can compare to Fubini-type theorems in HCN05I 
Thm 3.7] and [Mis08, Thm 1.13.1]. Our result extends these Fubini theorems to inte- 
grals of Hida distributions. 

Theorem 3.9 (Fubini-Tonelli theorem). Let 

F s ,r = ^2 c /3( s > r )^/3 
/3eA 

be an (S)* -valued process indexed by (s,r) £ R x [0,t]. If, for each ((3,k) pair, 
cp(s, r)Mnik{ s ) iJ bounded above or below by an L ([r, t] x [0, t\) function, then 

J J F s ^ r (uj) dB^ dr = J* QT F s , r (w) dr^j dB^ . (3.24) 

The equality in (13.241 ) is in the sense that if one side is in (S)*, then the other is as 
well, and they are equal. If in addition, 

(F s , r l [rA (s)oW H (s)) (s,r) e L\R x [0,t]), (3.25) 

then both sides are in (S)* and cg(s, •) £ L [0, s]for a.e. s S [0,i]. 

Proof. Unraveling the above definitions we have 

f [ F s , r {u)dB^ dr (3.26) 

Jo Jr 

F s>r l [rA (s)oW H (s)dsdr 



I I I 5Z^( s . r )l[r,t](s)^M ) o VM^(s)-H £ ( t) (w) dsdr 

Jo Jr \peA J \k>o J 

I I | ^ c p{ s ^ r ) l [r,t}{s)M H i k {s)li f}+£( k)(uj) J dsdr. 

Jo Jr y^gAjfeeN / 
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Denote the right-hand side above as G, an (5)*-valued integral. By Lemma [3T2l such 
integrals are characterized by their action on (5): 

ft 



((0,4-)) 



E 

/3eA,fceN 

E 

/3eA,fceN 



((c (s,r)l [r . t] (s)M H Ck(s)n^ +£(kh ^))dsdr (3.21) 

/3eA,fceN 
t 



((cp(s, r)l[ r>t ] (s)M H i k (s)Up +£ {k) ,il>))ds dr 
{{cp(s, r)l[ 0)t ] {s)M H ^k(s)'H /3+e (k) , ip)) dr ds 



■ /3eA,fceN 



,'0)} ds 



where equality in the second line follows since for any G A there are finitely many 
pairs (/3, A;) such that /3 + = j3' (recall also that the Hp are orthogonal). The third 
equality follows from Tonelli's Theorem for real-valued functions which is possible 
due to our hypothesis. If in addition (13.251 ) holds then both sides of this equality are in 
(5)* by (13 - 12b - The final equality follows from Lemma [ 
Teasing apart the right side of (I3.27t gives 



G 



E 

V/3SA 



cp(s, r)dr 



l m { S )Up{u) o J2 MH ^(s)n eik) (Lu) ds 



K k>0 



j* (j' F s , r (uj)dr) dBf 



as needed. 



(3.28) 
□ 



4 The Wiener chaos decomposition of DSLT 

In this secion we start by calculating the Wiener chaos expansion for a' t (0) defined in 
(1 1 - 10b - In the process, we obtain a new proof of the existence of a' t (0) for H < 2/3. 
Later in the section we will adapt the arguments used in obtaining the Wiener chaos 
to show existence of a' t (y) for all y G R. To reduce notation, in this sequel we write 
H = L 2 (R). 

Theorem 4.1. For H < 2/3, a' t (0) is in L 2 (V) and its Wiener chaos decomposition 
is 



oo 

4(0) = Y, I ^-i(9(2m-l,t)) 

m=l 



(4.1) 
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where g{2m - 1, t) G H 02 ™" 1 and 



s(2m-l,i) = g(2m - l,t;v 1 , . . . ,v 2m -i) 



.y\.L) 



(m-l)!2 m - 1 V2^7 Jo ( s - r )^( 2m - 1 )+ 1 

Proo/ Recall that Bp := (to, Mff lr )t ]) so that in particular, its Malliavin derivative 

is DBf = M H l[ ,t)- 

It is not hard to see that f' e (Bp - Bp) is in n fceN B fc ' 2 , thus by Stroock's formula, 
the n-th integrand in the chaos expansion of a' t £ (0) is given by 



1 ft 




o Jo 



(s - r) 2H - l V[D n f' e {BP - B?)} dr ds (4.3) 

-A \ (s-rf H -^[(-^f £ )(BP-BP)]YlM H l [riS] ( V] )drds. 
J J ■ i 



As in Section |21 we write 



Thus, 

tt „ , , £J 117,-, fc 



f e (x) = £- [ e^p n e-^l 2 dp. (4.4) 



'71 -1-1 

7 ! ' ' p n+l e -{(s-r) 2H +e)p 2 /2 dp 



2tt 



„-n+l 



27T((,s-r) 2 ^ + e)W 2 )+ 1 

i n+1 \/2^ (n + 1)! 

2vr((s - r) 2H + e)WW& 2(™+ 1 )/ 2 ((n + l)/2)! ' 



(4.5) 



if n + 1 is even, and if n + 1 is odd. Setting n = 2m — 1, it follows that the chaos 
expansion for ol t £ (0) is 

°° ( ( l) m /"* f s (s r )2i7— 1 2m— 1 \ 

((m-l^i^F/o /o FRT^ n M H l [r , s] ( Vj )drdsj . 

(4.6) 
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We now need to show that as e — > 0, the above converges in L 2 (V) to (14- 1 b - We 
will apply the following lemma adapted from [NV92|, which is a consequence of the 
Dominated Convergence Theorem. 

Lemma 4.2. Let F £ be a family of L 2 (V) random variables with chaos expansions 
F £ = X^^Lo In(fn)- lff or each n ' fn converges in H® n to f n as e — > 0, and if 

oo oo 

^supE[/ n (/^) 2 ] = £sup{n!| < oo, (4.7) 

n=0 £ n=0 6 

then F £ converges in L 2 (V) to F = Y^^=o ^n(fn) as e — > 0. 

We note that this argument has also been used in [HN05] and [Marl2|. To apply 
the lemma here, we calculate the L 2 {V) -norms of the chaos expansions and show 
they are bounded uniformly in e. Recall that T> t = {0 < r < s < t}. If we let 
g(2m — 1, t, e) be the integrand of hm-i in ( 14.61 ), we have 



E[/ 2m _i( 5 (2m- M,e)) 2 ] = (m - l)!||s(2m - 1, t, e)\ |^® 2m -i 



(2m- l)!((2m)!) 2 f (s - r) 2 ^ 1 ^' - r') 2H - 1 



27r[(2m - l)!] 2 (m!) 2 2 2m J V 2 (e + (s - r )2# )(2m+l)/2( e + jy _ r /)2H)(2m+l)/2 

2m- 1 \ 

J| (Mi J l [riS ]('yy),MH-l[ r ./ iS /](uj))Hrf'yj drdr'dsds'. 

(4.8) 



ira-l 



Maximizing by setting e = and using (A. 10) in [EVDH03] and the notation in (12.6I ). 
this is 



m(2m)! f (s - r) 2H "V - r') 2 ^V 2m_1 , , , 

/ " ■ ZTTo Lio dr dr ds ds ■ ( 4 - 9 ) 

J V 2 \m+l/2pm+l/2 



vr(m!) 2 2 2m J V 2 
Let 7 = /U 2 /(Ap). The L 2 (P)-norm of <@2) is then 



1 /" m(2m)! 7 m \ (s-r) 2 ^ 1 ^' - r') 2 ^ -1 drdr'dsds' 

However, by the generalized binomial theorem, 

7 ^ m(2m)\j m 



2(l-7) 3/2 ^ (^) 2 2 2 ^ 
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Thus, the L 2 ("P)-norm of (@~6]> is 

1 r 7 ( s _ r )2H-l( s /_ r /)2H-l 1 r _ r )2-ff-l( s / _ r /)2H-l 



27T^ f (l- 7 )3/2 MV ^p 2^ 7^2 (Ap-/,2)3/2 

By Lemmaim this is finite if # < 2/3. □ 

Remark: One might think of fBm as an isonormal Gaussian process W : H' — > 
L 2 {V) on the space H' = L 2 H {M) which may contain distributions. Then 

(l[0,t]> ![o,s])h' := (AfFl[0,t].^Hl[0,(])L 2 (E) (4-12) 
= i(^ + s 2if_| f _ s |2if ) 

Using this so call "twisted" inner product Hilbert space, the isonormal Gaussian pro- 
cess gives Bf = (u, lr 0) t] )tw so that DBf = ljo t] ■ Comparing this with the above, 
we see that the twisted inner product incorporates the operation of Mh in (oj, Mff/) 
into (oj, f)tw When / is a step function, essentially nothing but notation has changed, 
which can be verified by the use of the D M operator in Example 6.4 from [EVD H031 . 
One may then write 

rt j'S r g _ r \2H-l 2m - 1 rt rv f g _ r \2H-l 

Jo Jo (( S -rm^ + m g h,s ] (v j )drds = ^ ^ - _ r)2H){2m+1)/2 dr ds, 

(4.13) 

where 

V = Vi V ... V V2m-1, 

V = V\A...A V 2 m-1- 

It is then straightforward to verify that Proposition 14. II simplifies to the chaos expan- 
sion given in [Marl2] in the case H = 1/2. 

We now use the methods in the above proof to show L 2 (V) convergence for a' te (y) 
as e — > 0. 

Proposition 4.3. For H < 2/3 and any y £ R, a' t (y) is in L 2 {V). 

Proof. We may follow the proof of Theorem 14.11 except that in place of (14.51 ) we have 



m-^f*){Bf-B?-v)]\ = |^ jne ip{B? - B?) y py P n+l e- ep/2 d P 

1 

2 7r (( s -r) 2 ^ + e)K2)+i 



< 7; 77 2 , XTw^n / \ P \ n+1 e-^/*dp. (4.14) 
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We aim to apply Lemma l4~2l again. The arguments from Theorem 14. II show that the 
sum of the odd terms in (I4.7t converges. However, we can no longer argue that the even 
terms are 0, as we did before. Instead, we must use the identity J" R \p\ n+1 e~ p I 2 dp = 
2 n / 2 (n/2)!, valid for even n. 
Replacing (14.8I ). we then have 



E[I 2m {g(2m,t,e)) 2 } = (2m)\\\g{2m,t,e)\\ 2 2 

(2m)!(m!) 2 2 2m f (s - r) 2 ^ 1 ^' - r') 2 ^" 1 



(2m!) 2 J V 2 (e + (s- r) 2 ^)"^ 1 ^ + (s' - r') 2 ^) m+1 

„ 2m \ 

/ \\(MHh r ,s\{vj),M H l\r' ,s']{vj))ndvj dr dr' ds ds' . 
Mi™ H[ / 

(4.15) 

We proceed through steps (14.91 ) and (14.1 OK setting e = and 7 = p 2 /(Xp), in order 
to reach a bound on the even terms of 



(m!) 2 2 2m 7 m \ (s-r) 2 ^-V - r') 211 ' 1 dr dr' ds ds' 

*J*\Z-w-) v ■ <416) 

We now use the following identity and bound, valid for < 7 < 1 and with K a 
positive constant: 

» ( m! )2 2 2m 7 m ^ 7 ^/T^ + ^7 sin" 1 ( ^7) 

^ (2m)! ~ (1-7)3/2 

K 7 4 1?) 

- (1 _ 7 )3/2 

< ^ 
" (l- 7 )3/2- 

Inserting this and the expression for 7 into (14.101 ) gives a bound of 



if f p{s-r) 2H - l {s' - r') 2H - 1 dr dr' dsds' 



7T J V 2 (! _ g)3/2 A 3/2 /) 3/2 

K [ ^s-rf^js' -r') 211 - 1 
k L (Ap- M 2 ) 3 / 2 



(4.18) 



dr dr' ds ds' < 00. 



By LemmaH this is finite for H < 2 /3. □ 
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5 The Tanaka formula 

The following is what we have referred to as the Tanaka formula for the DSLT process. 
Theorem 5.1. ForO < H < 2/3, the following equality holds for all y and t in L 2 (V): 

Ha' t {y) + \sgn{y)t = I Lf-^dBf - \ f ' sgn(Bf - Bf - y) dr . (5.1) 
1 Jo 1 Jo 

Proof. Let f £ be defined as in Section [2] and let 

F £ (x) = / f £ {u)du= / f{u)du. (5.2) 
Jo Jo 

We apply Theorem 13.81 (Ito's formula) to F £ {x — y) using the fractional Brownian 
motion Bf — Bf, and then integrate with respect to r from to t to get 

f F £ (Bf - Bf -y)dr- tF £ (-y) = (5.3) 
Jo 

f f f £ (Bf - Bf - y) dBf dr + H f f f e (Bf - Bf - y)(s - r) 2 ^ 1 ds dr. 

Jo Jr Jo Jr 

Note that F £ (-y) -> -| sgn(y) as e -)• 0. Also, \F £ (-)\ < 1/2 for all e > 0, so by 
dominated convergence the integral term on the left side approaches 

~ f\gn(Bf-Bf-y)dr 
z Jo 

as e — > 0. We now want to apply Theorem I3T91 to the first term on the left side to get 

f f f £ (Bf - Bf - y) dr dBf . (5.4) 
Jo Jo 

To show (15.4I ). we use the following Wiener-I6 chaos expansion obtained from 
Stroock's formula: 

-E[(^_/ B )(Bf -Bf -y)](M*l M )®») • (5.5) 

As stated earlier, the Hermite chaos refines the nth Wiener-Ito chaos in terms of the 
Hermite orthonormal basis of H® n . In particular, the coefficient cp(s, r), with |/3| = 
n, in the Hermite chaos expansion of f e (Bf — Bf — y) is given by 

1 d n 

c p (s,r) = _E[(— / 6 )(£f - Bf - y)]{^, (M H l [r , s] )® n ) M ® n . (5.6) 
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Using (14.141 ) and the fact that M H £k is bounded (see Lemma 4.1 of (EVDHOl), one 
can now easily verify the conditions of Theorem |3.9| for f £ (B^ — — y) for e > 0. 

For fixed s, as e — ► the inner integral in (15.41 ) converges to Lf s y in (S)* 

by Proposition 10.1.13 in |BH0ZO8| . Now, if L? "~ y o W H {s) is integrable in {Sf , 
then by the arguments of Proposition 8.1 in [HKPS93], ( 15.41) converges in (S)* to 

Jo Lf s v oW H (s) ds. In other words, the equality in (15.11 ) is valid as long as one side 
or the other is in (5)*. However, by Proposition @3J a' t (y) G L 2 (V) , for H < 2/3. 
Thus for such H, <ED holds in L 2 (P). 

□ 

Remark: One-dimensional fBm has an L 2 (V) local time for any < H < 1 

(see MBH0ZO8B ), but it is not clear whether or not J* Lf s ~ y dBf is in L 2 (7 : ') for 
H > 2/3. As stated in the conjecture of Section|2j we suspect that it is not. Of course, 
a positive answer to the conjecture does not rule out the possibility that 

/ Lf"- y dB^ G (S)* for H> 2/3. (5.7) 
J o 

If (15.71) were indeed true, then the DSLT of fBm would also be well-defined in (5)* 
for all H G (0, 1). For H < 2/3, another open problem is to prove joint continuity, in 
y and t, of a' t (y) + sgn(y)i. One approach to proving this is to use the explicit chaos 
expansion for this integral (see Theorem 14 .11 ) combined with Definition 13 .7 1 



A Second moment bounds for fBm 

Recall 



A 


:= Var(5f - 




= | S -r| 2 ^ 


P 


:= Var(£^- 


- B$) 


= \s'-r'\ 2H 


n 


:= Cov(5f 




B^ - B^) 




= -Qs'-r 
2 VI 


_j_ 


s-r'\ 2H -\t 



(A.l) 



\2H 



r\ 2H ). 



The following useful bounds appear as Lemma 3.1 of [HuOl ]. As before, we let 
K denote a positive constant, which may change from line to line. 



Lemma A.l. (i) Suppose r < r' < s < s'. Let a 

Then 



r, b 



r', c 



s. 



\p-H 2 >K {(a + b) 2H c 2H + a 2H (b + c) 2H ) . 



(A.2) 
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( ii) Suppose r < r' < s' < s. Let a = r' — r, b = s' — r' , c = s — s' . Then 

Xp-p 2 > Kb 2H (a + b + c) 2H . (A3) 
( Hi) Suppose r < s < r' < s'. Let a = s — r, b = r' — s, c = s' — r'. Then 

\p-V 2 >K(a 2H c 2H ). (A.4) 

Unfortunately, (ii) is false, which can be seen by noting that if r' \ r, s' /* s 
then the left side of (1A.3I) approaches while the right side does not. However, (ii) 
can be replaced with the following, which suffices in every instance in which we have 
seen (ii) applied. 

( ii ') Suppose r < r' < s' < s. Let a = r' — r, b = s' — r' , c = s — s' . Then 

Xp-ti 2 >Kb 2H (a 2H + c 2H ). (A.5) 

Proof. We will follow the method of proof given for (ii) to arrive at (ii 1 ). We use the 
property of fBm known as local nondeterminism, which implies that the following is 
true for to < • • • < tj '■ 

Var ( J>;« - <j) > K^WiWU ~ ti-i? H ■ (A.6) 
\i=l / i=l 

In our case, this implies 

Var (u(£f - B?) + v(B$ - B$ )) (A.7) 
= Var [u(B$ - B* ) + (u + v)(B$ - B$) + u{B? - B$)) 
> K (u 2 (a 2H + c 2H ) + (u + v) 2 (b 2H )) . 

However, we also have 

Var (u(Bf - B? ) + v(B$ - Bf, )) = u 2 \ + 2uvp + v 2 p. (A.8) 
Combining (IA.7I ) and (IA.8b yields 

u 2 (\ - Ka 2H - Kb 2H - Kc 2H ) + 2uv(p - Kb 2H ) + v 2 (p - Kb 2H ) > 0. (A.9) 
The discriminant of (|A.9b therefore satisfies 
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4(/x - Kb 2H f - 4(A - Ka 2H - Kb 2H - Kc 2H )(p - Kb 2H ) < 
which gives 

Xp-p 2 (A. 10) 

> K{(a 2H + b 2H + c 2H )p + b 2H X-2pb 2H ) 

—K 2 (b 2H (a 2H + b 2H + c 2H ) - b m ) . 

By the Cauchy-Schwarz inequality, 

p<Vx~p<^. 

Using this, together with A = (a + b + c) 2H and p = b 2H allows us to reduce (IA.10I ) 
to 

\p-p 2 >(K-K 2 )b 2H {a 2H + c 2H ). 
The result follows by replacing (K — K 2 ) with K. □ 

Proof of Lemma 122} Recall that we must show 

f u ( s _ r )2H-l( / _ M2if-1 

/ drdr'dsds' <oo, (A.ll) 

for < il < 2/3. We will split the range of integration into the three regions de- 
scribed in in Lemma lATTl with a, 6, c defined accordingly on each region. 

Case 1: r < r' < s < s'. Using s — r = a + b, s' — r' = b + c, and (IA.2b we see 
that the contribution of this region to (lA.llt can be bounded by 



J m > ((a + b) 2H c 2H + a 2H (b + c) 2H fl 2 

f u(a + b) 2H ~ l (b + c) 2H ~ l 

~ K J [0>t]3 ( a + 6 )3 g /2 c 3 g / 2a3g /2 (b + c)3g /2 d «* dC < A " 12 > 



/x da db dc 

[0 , t] 3 (a + b) l - H / 2 c 3H / 2 a 3H / 2 (b + c)i-^/2 ' 



Using an idea which appears in MHuOl] , we write 
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2^ = {a + b + c) 2H + b 2H -a 2H -c 2H 

= 2H(b + c) / (a + ib + c^-'du + b 211 -c 2H . 
Jo 

We can bound the integral in (IA. 131 > by replacing the integrand by its maximal value. 
This gives a bound of \fi\ <K((b + c)a 2H ~ l + b 2H + c 2H ) for H < 1 /2, and a bound 
of |//| < if((6 + c) + b 2H + c 2H ) wheni? > 1/2. 
If H < 1 /2, we bound ([021 by 

^ /' ((6 + c ) a 2g-i + 6 2 g + c 2 g) 

y [0it] a ( a + 6)1-^/2^/2^/2^ + c) l-ff/2 dad6dc - 

We need only show that each of these three terms is integrable, and we may obtain 
bounds by replacing nonnegative powers of (a + b) in the denominator by powers of 
either a or b, and likewise for (b + c). When H < 1/2 we will also use (a + 6) 1- ^/ 2 > 
a i-2H b 3H/2 i n tegrability follows from: 



(b + c)a 2H - 1 K K 

( a + b y-H/2 c 3H/2 a 3H/2( b + c y-H/2 ~ ( Q + b y-H/2<3H/2 a l-H/2 ~ b l-H/2f3H/2 a l-H/2 ' 

b^_ . Kb 2H _ K 

( a + b y-H/2 c 3H/2 a 3H/2( b + c y-H/2 ~ b 3H c l-H/2 fl l-H/2 ~ & H c l-H/2 a l-H/2 ' 

a; < k 

^ a + b y-H/2 c 3H/2 a 3H/2^ b + c y-H/2 ~ ( a + b y~H/2 c l-H a 3H/2 ~ b l-H/2 c l-H a 3H/2 ' 



In each case we obtain an integrable function of a, b, c, i.e. the exponent of each 
variable in the denominator is less than one. 
For 1/2 < H < 2/3, we have 

(b + c) K K 



^ a + b y-H/2 c 3H/2 a 3H/2 {b + c y-H/2 ~ ( fl + b y-H/2 c 3H/2 a 3H/2 ~ b l-H/2 c 3H/2 a 3H/2- 



b 2H R 

< 



( a _|_ b y-H/2 c 3H/2 a 3H/2( b _|_ c y~H/2 — b 2-3H c 3H/2 Q 3H/2 ' 

while the third term follows exactly as in the case H < 1/2 (note the second inequality 
in dA. 15b used a bound which is not valid when H > 1/2). This completes Case 1. 

Case 2: r < r' < s' < s. Following IH uOll we can write 
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o ( i u\2H 2H . /. , \2H 2H 

L\x = (a + b) —a + (b + c) — c 

= 2Hb f \{a + bu) 2H - 1 + {c + buf H - 1 )du. <A ' 16) 

JO 

Replacing the integrand with its maximum over the region gives us a bound of n < Kb 
for H > 1/2 and fi < Kb{a 2H ~ l + a 211 ' 1 ) for H < 1/2. 

First consider if < 1/2. In this case, s — r = a + b + c and s' — r' = b. Also note 

that 

(a + c) 2ff < (2ma X (a,c)) 2H = 2 2ff ma X (a,c) 2/f < 2 2ff (a 2ff + c 2ff ), 

thus (a + c) 2H and (a 2H + c 2H ) are equivalent up to a constant. We bound the contri- 
bution of this region to (IA. lib , using (IA.5I ). by 



, 6(a 2 ^" 1 + c 2H - l ){a + b + c) 2H - l b 2H - 1 , ^ J 

ii / -^ttt -^Tj da db dc 

'[o,tp b 3H (a + c) 3H 

< A / -JT- — ■ — -T777 aa ao ac 



[0,t]3 



,2H-iy2H-l 



<K 7777 — "—^dadbdc (A.17) 

< K da dc 

y [0it] 2 (a + c) 3 ^ 



'[o,t] 

For 1/2 < < 2/3, we have 



< K f a 2H - l {l + -^r) da < oo. 



/ 

■/[0,t] 



3 b 3H {a + c) 3H 

r ( a -i-h-i- c )2H-i 
= K Z dadbdc (A.18) 

< K f dadbdc 



This completes Case 2. 

Case 3: r < s < r' < s'. Following IH uOlU we write 
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2fl = ( a + b + c) 2H + b 2H - (a + b) 2H - (b + c) 2H 

= 2H(2H-l)acf [ (b + vc + ua) 2H ~ 2 dudv. 
Jo Jo 

By Young's inequality, applied twice, for a, (3 > with a + B = 1 we have 



(A. 19) 



(6 + + ua) > Kb a (vc + uaf > Kb a {vcf/ 2 {uaf' 2 . (A.20) 
This combines with (1A.19I ) to give 

|M| <K(ac)^ H -^ +l b 2ai > H - l \ (A.21) 

Using s — r = a, s' — r' = c, dA.211 ), and (|A.4b shows that the contribution of this 
region to (IA. lib is bounded by 

f (ac)^ H -^ +l b 2 ^ H - l )a 2H - l c 2H - 1 „ , 
K I „„ „„ da db dc 

m (A.22) 

da db dc 



K 



W b 2a i l - H ){acY+ H ^-P)' 
Note that 1 > 2H{1 - H) so that w^H) > H - We 

may therefore choose a such 
that H < a < 2 (i-H) implying P < 1 — -ff. The exponents in the final expression 
in dA.221 ) can therefore be bounded by 2a(l — H) < 2( ^-1 )(1 - H) = 1 and 
/3 + - /3) < (1 - F) + H = 1. We conclude that dX22l is finite, which 
completes the final case for the proof of Lemma [2721 □ 
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